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1 Introduction 

The Evans function is a tool for assessing the stabihty of travelhng waves solu- 
tions for partial differential equations. 

A recent paper [2] analyzes the order reduction experienced when evaluating 
the Evans function numerically. The details of some lengthy calculations were 
excluded from that paper for clarity. The purpose of this technical report is to 
make these details publicly available. This report is not intended to be read on 
its own; the reader is referred to [2^ for background and references. 



2 Setting 

We consider scalar reaction-diffusion equations of the form 

Ut = Uxx + /(")■ 

Let u{x, t) — u{S) with ^ = x — ct be a travelling wave solution of this equation. 
A linear stability analysis of this travelling wave leads to the eigenvalue problem 



where 



The limits of A as ^ 



= ^(^;A)2/, 

1 ■ 

A-/'(u(e)) -c 
±oo are given by 



A(?;A) 



(la) 

(lb) 



1 

A-/'(ii±) -c 
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Furthermore, the eigenvalues of (A) are 

A^™ = ^ (-C + v/c2+4(A-/'(fi_))) and 
- ^ (-C - v/c2+4(A-/'(fi_))) , 
where u- = hm^^_oo ii-iS,)- Similarly, the eigenvalues of A^{X) are 
= i (-C + v/c2+4(A-/'(fi+))) and 
A^?' - ^ (-c-v/c2+4(A-/'(fi+))) , 

where u+ = lim^_^+oo The corresponding eigenvectors are and 

respectively. 
To define the Evans function, assume that 



ReA>max(/'(u_),/'(u+)) 



ImA 



(2) 



This condition ensures that and fj}^^ have positive real part, while fJ3 and 

[21 

/xL have negative real part. 

The differential equation ([1]) is linear, and hence its solutions form a linear 
space. Let y- be the solution which satisfies 



as ^ 



(3) 



Condition ^ implies that this defines y- uniquely, that satisfies the bound- 
ary condition — s- as ^ ^ —oo, and that any solutions satisfying this 
boundary condition is a multiple of y_. 

Similarly, we define as the solution satisfying 



y+iO-e^v{^^^lk) 



1 

[2] 



as ^ — > oo. 



(4) 



The Evans function is then the function D defined by 

i^(A)=det[y_(0)|y+(0)]. 
We are interested in computing this function. 



3 Asymptotics near infinity 

In this section, we study the behaviour of 13(A) as 
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3.1 The solution satisfying the left boundary condition 

Define a transformation y_ y_ by 



where 



1 



1 

,[2] 



and _B_ = 



exp{^l^lk)B-y-{0, (5) 
(6) 



The differential equation ^ transforms to 



at; K_ K_ 



1 



1 



where 



(7) 

(8) 
(9) 



V-iO = fiHO) - /'("-) and K_ = Vc2+4(A-/'(fi_)), 
and the boundary condition ^ becomes 

u-{^) 1 and w-(^) ^0 as ^ — > —oo. 
Now, suppose that u_ and v- can be expanded in inverse powers of k_: 

We now substitute these expansions in ([7]) and equate the coefficients of the 
powers of 

• At 0{k-), we get = Vq , so Vq is identicaUy zero. 



At we get (u^ )' = and (v^)' 



The first equation, together 



with the boundary condition ([9|), implies that Uq =1. It follows from the 
second equation that is identically zero. 

• At 0{kZ^), we get 

{ui Y = -ip-iO{uo +Vo ) and {v^)' = tp^{^) {u^ + v^) - ■ 

Substituting what we found before yields 

(mJ")' = -V9_(0 and = (/3_(O-w^. 

Hence, u^{£,) = — <f>-(^) and V2 {£,) = </'-(0; where 



(x) dx. 
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• At ©(kI^), we get 

{^2)' = -V-iOiui +vT) and (u^)' = (^_(^) (uj^ + ) - 5;^. 

Substituting and in the first equation yields {112)' — f-iO ^-iO-i 
and hence, 



U2{0 



(p-ix) / (p^{y)dydx 



— oc J —00 



Summarizing, we have the foUowing series expansions for the transformed soki 
tion: 



(10) 



3.2 The solution satisfying the right boundary condition 

The computation for the solution y+ satisfying the right boundary condition ^ 
is analogous. Define the transformation y^ — > y+ by 



y+{0=^Ml^^lk)(u+{0 
where 



= cxp(A*V''Os+y+(e) 



y+ 



u+ 



and B 



+ 



The differential equation ([T]) transforms to 



(11) 



where 



^+(0 - /'("(O) - /'(^+) and = Vc2+4(A-/'(7i+)), (12) 
and the boundary condition ([4]) becomes 

u+{£,) and w+(C) ^1 as ^ ^ 00. 
Expand u+ and 7}+ in inverse powers of k+: 

i2+(C; AJ+) - 4(0 + '^;'4(e) + >^+'u+io + ^+'u+{o + o{k-^'), 
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Substituting these expansions in the transformed differential equation PT|) and 
equating the coefficients of the powers of k+ yields the following equations: 



= 4, 

{u+y^-^+iO {ut + v+) + ut, 



= 0, 
(^o+)' = 0. 



Solving these equations in the same way as in the previous case yields the 
following series expansions for the transformed solution: 



(13) 



where 



$+(^) 



(p-)-(a;) da;. 



3.3 Asymptotics for the Evans function 

The Evans function is obtained by evaluating both j/_ and ?/+ at ^ = and 
taking the wedge product (which equals the determinant of the 2-by-2 matrix 
having these vectors as its columns). This yields 



D{\) = y^(0) A 2/+(0) = (S-y-(O)) A (B+y+(0)) 

1 1 1 rw„(o)i r 1 

_i(K_ - c) -i(K_ + c)\ [w_(0)J - c) 

+ 1 (^_ + n+){v_{Q)u+{Q) - u_{Q)v+{Q)) 
It follows from (jTU]) and that, as |A| oo, 



1 

-i(K+ +C) 



u+(oy 

v+iO) 



V- 



u+ 



OK') 



O(A-i), 



O(A-i), 5+ = 0(l). 

Furthermore, we have k_ — K4- = 0{X-^/^). Hence, 

D{X) = -i(K_ + «;+)w(0)w+(0) + O(A-i). 



Finally, 
where 



-2A1/2 + $ - iA-i/^l^*^ - 2/'(it_) - 2f'iu+) + c') + O(A-i), 



$ = $_(0) + $+(0) 



(/3_(a;) dx ■ 



^+{x) dx. 



5 



4 The exponential midpoint rule 

The exponential midpoint rule (or second-order Magnus method) for solving the 
differential equation y' — A{^)y is 

yu+i = exp {hA{£^k + \h))yk. 

We assume that the step size h is fixed, so that ^fe = + ^h. 

The transformation ^ changes the recursion for the exponential midpoint 
rule to 

yk+i = exp(-^'^'^„+i) BZ^ yu+i 

= exp(-^We/c+i) BZ^ exp {hA{£,k + jh))yk 
= exp(-^[^'^fe+i) BZ^ exp [hA{^k + \h)) exp(^'^'^fe) B_ yu 
= exp(-V™) SZ^ exp (/iBzMfe + \h)B^) eMt^^-kk) Vk 
= exp + \h))yk, 



where 

Using (jlbp . we find that 



A.iO^BZ'AiOB^-fi^lh. 



with if- and k_ as defined in 

The transformed recursion for the exponential midpoint rule is the same as 
the exponential midpoint rule applied to the transformed equation y' — ^(^) y, 
cf. ([7]). The reason for this is that the Magnus method is equivariant under 
transformations such as (O. 

4.1 The local error 

The local error for the exponential midpoint rule is defined by 

L- = exp {hA{^k + ^h)) y{^k) - 2/(6+i), 

or, in transformed coordinates, 

L- = exp {hA^i^k + ^h)) - 

We compute the matrix exponential of hA by diagonalization. The eigenvalues 
of are 



6 



and 



\2 = -^h(^K^ + ^nl-A^^iO) 



The corresponding eigenvector matrix is 

V 

with 



1 1 

Vi V2 



- k^/k^ ^Aifi) -1, 



2ip 



where we are writing k for k_ and (p for ip^{^). Its inverse is 



V2 — Vi 



V2 



1 



We have hA = V[^^^ a° ]^"^ and thus exp(/iA) = 1/['=''p^^ cxpAs]^"^- How- 
ever, A2 ^ —hn so that exp A2 is exponentially small as |k| — > 00 if k is restricted 
to lie in a sector of the form | arg^j < ^tt — e. Under this assumption, 



exp(/i^) = V 

exp Ai 



exp Ai 
e.s.t. 



V2 
V1V2 



-1 

-Vi 



e.s.t. 



V2 — Vl 

exp(— ~ 4(p)) 



1 — ^ 

K 



hip h'^ip^ 



2/^2 



-3^ 



e.s.t. 



(14) 



where e.s.t. stands for exponentially small terms. 
Hence, using the exact solution pI7|) . we find that 

exp(/iA(a + i/i))y(6) 



1 - 



2k? 
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We thus arrive at the following expression for the local error: 
L- = exp {hA{Ck + ^h)) y{^k) - vi^k+i) 



1^ 



V{ik + \h)-ip(^k + i) 



-U~^hip'i^k + kh) + 0{K-^h^, K-^h) 



We need to assume that \k\ ^ h ^ for the last equality. 



4.2 The global error 

We write the local error as 



where 



Ik 



LI 



ik+h 



if- (x) dx - V- (6 + ^h)^ Oih"), 

5k = ip-i^k + hh) - V-{ik + h) = 0{h). 
The global error satisfies the recurrence relation 

Ek+i = cxp(/iA(a + \h)) + LI, = 0, 

or, in transformed coordinates, 



^^+1 = exp(/iA_(a + \h)) E- + L 



k ' 



E- 



0. 



The solution of this recursion is 



H-^Sk-l + 0{K-^h) 



(15) 



(16) 



This can easily be proved by induction. The case fc = 1 is trivial. Assuming 
that holds for a particular value of fc, we have, using (|14p and p^ . 



4+1 = exp(/iA_(a + \h)) 4 + LI 



1 + 0{k-^) 



K-'5k-l+0{K-'h) 
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which concludes the induction. 

Substituting 7^ and Sk back in p6)) . wc find that 



E7 



K"^ ((^_ (Cfc + i/i) - ^- (a + h)) + 0{K-^h) 
-1 ip{x) dx - hY^tl ^(e, + \h:)) + 0(«:-2/i2y 

where (^(^) = f'{u{0) differs from (p-{£,) by a constant. 
4.3 The solution on [0, 00) 

To compute tlie solution y+ satisfying the right boundary condition, we run the 
exponential midpoint rule backwards: 

yk+i = exp{-hA+{£,k - jh)) yk, 



where ^fc = S^a — kh and 



(17) 



with and as defined in (|12p . A similar computation as before yields 

„2 



exp(— 



and 

exp(/il+(efe-i/i))y(a) 



^+(...) 



where (/?+(. . .) stands for (p+{£.k — ^h). When we use this to determine the local 
error, we find 
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A standard induction argument shows that the global error is 



i^fi^i^k + ^h) - ifi^k + h)) + 0{K^^h) 



1 I f?o 



■'j 



4.4 The error in the Evans function 

The numerically computed value for the Evans function is the wedge product 
of the numerical solutions: 

^„um(A) - hvt^ (2/-(0) + E^) A (y+(0) + 

with k chosen such that = 0. Hence the error in the Evans function is 

Ed{X) ^ D{X) - D (A) 

= y_ (0) A E+ + E- A y+ (0) + E' A E+ 

= B-y-{0) A B+E+ + B-E^ A B+y+{0) + B-E^ A B+E+. 

Substituting _B_ and _B+, we find 

Ed = ^{n- - >^+){v-{0) [Eth - ^i-(O) [E+h + [E^hv+{0) 

- ^I+(0) + [E^h [E^h - [E^h [Eth 



+ + K+){v-{0) [E^]i - "-(0) [E+h + [E^h S+(0) 



(18) 



- [E-],v+io) + [E-h mi - [E^]i mh 

We now estimate all the terms in this expression and drop the ones of lower 
order: 

Ed = -~^{k^+k+){u^{0) [E+h + [Ek]iv+{0) + [E^h [E^h) +0{X-'h) 
= h V{]h+\h)~ / ip{x)dx + 0{\-^'^h'^), 

„■ M J — L 



■j = -N 



assuming that the differential equations are solved on the intervals [— L,0] 
and [0, L] with L = Nh. The final step is to apply the Euler-MacLaurin sum- 
mation formula (see e.g. [T]), which states that 



3=0 



f{x)dx+^h{f{0) + finh)) 



„o L2m+3 



(2m + 2)! 



-/ 



(2m+2) 



(0 
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for some ^ G [0, nh], where Bk denote the Bernoulli numbers. This yields 
Ed= ip{x)dx+ / (p{x)dx + ^h{ip{-L) + ip{L)) 

J-L JL-^h 

i=l ''^^>- 

Now, ip{^) decays exponentially fast to zero as |^| — > oo. So if we assume that 
L is sufficiently large, we can ignore all terms in this equation but the last one. 
We thus arrive at the final result, which is that the error in the Evans function 
is of order A~^/^/i^. 

5 The fourth-order Magnus method 

We repeat the computation in the previous section for the fourth-order Magnus 
method. This method is given by 

where [ • , • ] denotes the matrix commutator defined by [X, Y] = XY — YX and 
^l, are the Gauss-Legendre points 

Cfe = a + (5 - 1^3)/^ and = + {I + lV3)h. (19) 

After the transformation ((5|), the method reads 

yk+i = exp(ilfc)yfc 

with 

= ^h{A^iQ + A,ie,)) - 4h^[A.ia),A.{Q] 



= h 



ak 



























where ak and f3k are given by 

ak = U^-i^l) + ^-i^l)) and p, = -^h{^^iek)-^-iek))- (20) 

5.1 The local error 

The eigenvalues of Ofc are 

Ai = -^h(K^ - xk) = -h^nZ'xk + >^-\l + 0(«"')) 
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and 



A2 = -^h{K^ + xk) = -h(n^ - nJxk - kJxI + 0{kJ)'^ , 
where Xfc = "^fe ~ f^l 

Xk - ,J^l - 4(afc - Pi). 
The corresponding eigenvector matrix is 



V = 



1 1 

Vi V2 



with 



Vi 



2ak - + K-Xk , 2afc - - n-Xk 



Its inverse is 



2(K_/3fe + ak) 



and V2 



2(K_/3fe + afc) 



1 



1^2 -1 
-VI 1 



W2 - t^l 

As in the previous section, where we were considering the exponential midpoint 
rule, we have Ai —hn so that expAi is exponentially small (under the same 
assumption as before). Hence, 



exp(ilfe) = V 

exp Ai 



exp Ai 
e.s.t. 



V2 -1 
V1V2 -Vi 



V2 - Vi 

exp(-i/i(K_ - Xk)) 



+ e.s.t. 

ak - K^Pk 



li-Xk 



2ak - Kt 



li^Xk 



2{K^Pk - ak) 
K-Pk + ak 



K^Pk ~ ak 



exp(-^/t(K_ - Xk)) 

l^-Xk 



2ak - + K-Xfc 
2(K_/3fe - ak) . 



+ e.s.t. 



^{n'i + K^Xk) ~ ak K-Pk-ak 



K-Pk + ak 



■^{nt - K^Xk) - ak 



\_hXk. h'xl-2Pl 



+ 0{t,-J) ^ + 0{k-') 



Pk 



+ OiKZ^) 



Pk 



+ OiKZ^) 



(21) 



where e.s.t. stands for exponentially small terms. 
Hence, using the exact solution (fTO|) . we find that 



exp(flfe)y(^fe) 



^ + OiKZ^) 

K- 
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We thus arrive at the following expression for the local error: 



Lk = exp(f^fe) y{ik) - y{£,k+i) 



where 



lk= I f- {x) dx - h{ak - Pk) 

Uk 



and 



Pk = ^h^v'{^k + hh) + 0{h^). 



5.2 The global error 

The global error satisfies the recurrence relation 

Sfe+i = exp(nfc) E- + L- , E- = 0. 
The solution of this recursion is 

-I'l 



E7 



KZ'Pk-i + 0{KZ'h) 



This can easily be proved by induction. The key step in the proof is the following 
computation: 



^fc+i = exp(nfe^) 4- + 

'i + 0{kZ^) 0{kZ^) 
0{kZ^) OinZ^) 

'nZ^lk + 0{KZ^h^) 

KZ^Pk+0{KZ^h) 



KZ'Pk-l+0{KZ^h) 



+ 



>^--Y.Ul3+0{KZ^h^) 

KZ^Pk + 0{KZ^h) 
Substituting 7^ back in the formula for E^ , we find that 

{g'^-{x)dx-hj:'zl{c^k-Pi))+o{K-'h^) 

K-^Pk + 0{KZ^h) 

We now approximate 

0k = -^h{^.{ek) - ^-iek)) = ^h'<p'{^k - + o{h% 



E: 
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Thus, 



and therefore, 



fc-i 

^{x) + ^,h\^'{x)r dx - ^hj2i^{a) + 

^KZ^h''^'{ik-\h) + 0{KZ^h^,KZ^h) 



5.3 The solution on [0, oo) 

To compute the solution i/+ satisfying the right boundary condition, we run the 
same method backwards: 

Vk+i = exp(f^^) Vk 

with 



0+ = -\h{A+{ek)+A+{ek)) - 



1)1 



where 



6=eo-fc/i, ek=ik-{\-\^i)h and = - (i + iV3)/i. (22) 



We can write the matrix as 



ctk a "fe' 
K+ (3k ^ 



K+ 



where 

ak^U^+m + f+iek)) and /J^ = - ^+(^2^). 
A similar computation as before yields 



exp(— fi+) 



lK+ 



+ 0{nf) 1 



Pk_ 

K+ 



hiak - /3|) ^ h\ak-Pir-2Pl ^ 



K+ 



2k\ 



and 



exp(-0+)) y(^fe) = 



c&+UA.) + /Kn,-/j^) 
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When we use this to determine the local error, we find 



1 



ip+{x) dx - h{ak -01)]+ 0{k^%^) 



As in the previous section, we conclude that the global error is given by 



fc-i 



v{x) + ^^h\^'{x)f Ax - \hY,{^{eu) + fm 

3=0 



5.4 The error in the Evans function 



Substituting the results for the global error in p8)) , we find that the error in 
evaluating the Evans function is 



N-l 



j=-N 



As with the exponential midpoint rule, the Euler-MacLaurin summation for- 
mula can be applied to show that the term 

\h J2 (^(J^ + (I - l^)'') + Hj'^ + (i + - / ^i^) dx (23) 

j=-N •'-^ 

is negligible if L is sufficiently large. So, we find that the error in the Evans 
function is given by 



Er 



jL-h'^ {^'{x)rdx + 0{X-'/^h^). 



6 The fourth-order Gauss— Legendre method 

The two-stage Gauss-Legendre method for solving the equation y' — A{^) y is 
given by 

si = Aiek){yk + \hs, + {\-f)hs2), 
s2^Aiek){yk + i^ + ^)hs, + \hs2), 
Vk+i =yk + ^h{si + S2), 
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where and ^| are the Gauss-Legendre points, given in (fTO|) . As usual, we 
transform this to 

yfe+1 = yfe + \h{si + S2). 
6.1 The local error 

Substituting A_ and = y{^k) in the above formula and rearranging yields 

Vi\ _ , , /icri(/3i_ , haiipi _ (pi _ ipi _ 

1 + -; — sii + -; — S12 H S21 H S22 = Vki yk2, 

4k / k k k k 

hipi _ f hK h(pi\ _ haiipi _ / \ _ 



— yki - [K 

K \ K 



Vk2, 



ha2<f2- , ha2f2- , /, , hip2\ - , h(p2 _ _ (^2 _ <y52 



-sii H S12 + -L + -; — S21 + - — S22 = Vki = yki yk2, 



h(J2f2 _ , , / V2 \ - h(f2 _ 

Sll + /lfT2 K S12 -. S21 



hip2 
4k 


S22 


= Vki 


V2 _ 

= Vki 

K 


< 


T 


+ 1 - 


h(p2\ _ 






"^2 _ 

— Vki 

K 





yk2, 



where 

o-i = i-^N/3, cr2 = 3 + iV3, (^1 = (^-.(^^.) and <^2 = V'-(Cfc)- 
We assume that the unknowns Sij can be expanded in powers of k like 

K K 

Substitute this in the set of four equations above and collect like powers of k. 
• At order k, we find 



/is?2 + haiS22 = -yk2, 

h(T2sl2 + hs22 = -27fc2- 



Solving these equations yields 



_o 2^/3 _ J -0 2\/3 _ 

S12 = T— yk2 and S22 = —r- yk2- 
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At order we find 

- n 

jhsl^ + + haisl2 = 0, 

ha2sl2 + \hs\2 + sSJa = 0- 
Substituting and S22 and solving th.e resulting set of equations yields 

-0 -0 n -1 12(V3-1) _ _i 12(V3 + 1) _ 
Sii = S21 = 0> S12 = 2/^2 and S22 = yk2- 

At order k"^, we find 

+ j/«/?is?i + |:/iV5is?2 + haiipi^2i + /icri(/3iS22 = -ViVki - 'Piyk2, 
-jhipi^i + \hsl2 + s\2 - lhipi^2 - haxip\s%x + haxs\2 - hox^PxA.-! 

= ^iVki + v>iyk2, 

ha2<fi2^1 + /l<T2<P2S?2 + S2I + jhip2S^i + \h(fi2^2 = -^2ykl - V2yk2-, 

-hai(p2^i + ha2s\2 - ha2(p2S^2 - \hifi2^i + i/is?2 + ^12 - \h<fi2S^2 

= ^iVki + ^iyk2- 

We substitute all the known quantities in these equation. The values of 
and S21 can then be found from the first and third equation, respectively: 

s\i = -fiVki and s^^ = -(p2yki- 

The second and fourth equation become 

jhsl2 + h(Tisl2 = fiVki - 12(\/3 - l)h~'^yk2, 
h<T2s\2 + \hsl2 = 'P2yki + 12(V3+ l)h^^yk2- 

The solution of this system is 

s?2 = ^(3^1 + (2\/3- 3)(p2)yfci - ^(2^3- 3)yfc2, 

^2 = -JliC^^ + 3)^1 - ^V2)yki + ^(2\/3 + 3)yfe2. 

At order the first and third equations are 

+ \hipisl^ + \h(pisl2 + haiipis\^ + h(7iLpis\2 = 0, 
h(72ip2s\i + ha2V2s\2 + .S21 + \hip2s\x + \h'p2s\2 = 0- 
The values of s\i and s\i follow immediately: 

Sii = ^(iVi + (i - \^)^i^2)yki - '^^vm2, 

S21 = ^((i + gV3)v5i</52 + I'filjyki + ^^¥'2yfe2. 
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Collecting the results, we find that the stage values for the Gauss-Legendre 
method are 



Sll 



Sl2 



S21 



S22 



'Pi. , /i(3^? + (3-2V3)^i^2)_ 273(^1 _ , 
yfci H Vki j^^yk2 + 0{k ), 

2V3_ 12(\/3-l)_ 

— T—yk2 -\ r5 yk2 

, 3^1 + (2V3 - 3)(^2 _ 24(2V3-3) _ .3 



LP2_^ , /l((3 + 2V3)(pi</?2 +3V5^)^_^ ^ 2V3(y92^^ , /o/..-3^ 



— yfei 



-yki 



-yk2 + 0{K-'), 



2V3_ 12(V3 + 1). 



-yk2 



h^K 



-yk2 



(2\/3 + 3)(^i -3(^2_ 24(2^3 + 3) 



-yki 



yk2 + 0{K-'^). 



The result of doing one step is therefore 
yk+1,1 = yki + + S21) 



yki 



h{(fi+ip2)- , ^^(¥'1 + ¥'2)^ _ , \/3((/32-Vl) 



2k; 



4k2 



and 

yfc+1,2 = yfc2 + ^^•(Sl2 + S22) 



12 _ 

yk2 - -7—yk2 

tlK 



V3((^2 -fi)- , 72 _ _3 



We can write this as yk+i = ^1. J/fc with 



1 



2„,2 



K 2k2 

12/?fc 



12/3fc 

12 72 



where ak and /3/c are defined by ([20]) . It follows that 



1 - 



K 2k2 
/i(p_(a) + 12/3fe 



/ik2 



Finally, the local error is given by 

Lk^^kvi^k)-mk+i)^ 
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where 



ja, 



$_(^fc+i) - $-(Cfe) - hak 

/ ^{x)dx-\h{^{ek)+^{ek)) 



and 



and 



</^-(«) + — T </5-Ufc+l) 



6.2 The global error 

The global error is given by the recursion 
The solution of this recursion relation is 

fc-l k-l / j-1 



0. 



~ j=0 ^ i=0 



- j=0 

Indeed, assuming that the result holds for some value of k, we have 
^k+i=%E^+Lk 



+ 



^ E -=0 + E .toi^r - eu 4"'-) + oit^z'h^) 
«='E-=o^r+^(«-'^') 

and the formula for the global error follows by induction. 



19 



6.3 The error on [0, oo) 

The solution on the interval [0,cxd) can be computed by running the Gauss- 
Legendre method backwards: 

Vk+i = j/fc - ^^(si + S2), 

with as given in (fT7|) and £,k, ^1 ^^'^ Cfc ^ given in ([22]) . The global error 
can be computed as before, but here we will take a short-cut. If we comparing 
the matrix v4_ with A^ and the exact solution on the interval (—00, 0] with 
the exact solution on [0, 00), we find that they can be related by swapping the 
components 1 and 2, replacing ^ by — ^, and replacing the — subscript with 
a + subscript. Hence, the global error of the Gauss-Legendre method run 
backwards is 

+ j=0 + j=0 ^ i=0 ^ 

t ^(x)dx-i/i((p(ci) + (p(42)) 

6.4 The error in the Evans function 

The error in the Evans function is given by (fT8|) : 

Ed - - «:+)(c-(0) [E+]2 ~ S_(0) [E+], + [E^hv+{0) 

- [E-], u+iO) + [E-]2 [E+h - [E~]^ [E+]A 

- ^+(0) + [4"]2 \ElV - \ElV \Et\^ ■ 
Estimating all the terms, we find that 

Ed = -^{n- + «+) (m-(0) [E+h + [E^]i v+iO)) + 0{\-^'^h\ 



where 
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Let X denote the expression between the big parentheses. We need to evaluate 
this expression: 



^ \ + j=0 + j=0 ^ i=0 ' 

+ E £?- + ^ E (^r - 1: )) (i - 



$+(0) 



,4n 



1 



3=0 

N^l / 3-1 j-1 



j=0 \ i=0 i=0 

The sum + ^) same as expression (1^ . which appeared in 

the fourth-order Magnus method. As we discussed there, this expression is 
neghgible. Substituting the values of and ij'^, we find that 

N-l / J-1 
j=0 \ i=0 



i=0 



Exchanging the double sums yields 
^ A 



X E ^7" ( 5^7" + *-(-^ + + 'J'+w + ^ E "7 ) 

\j=0 ^ i=3 ^ 

+ ^7^ (5^7^ + (i - + (0) + ^ «+ n . 

Now, L°j'~ was defined as 

= + j/i + /i) - $_(-L + j7i) - /la", 

and thus we have 

N-l N-l 

hY. ai^ *-(0) - '^-{-L + Jh) + ^ Z;-. 
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Using this expression, and its equivalent for ^ • a^ , we find that 



= - A E (^-(0) + ^+(0) + + E 

+ 1]'+ u (0) + *+(o) + + x; 4"'+) ) • 

We know that Z"'^ = Oih^), which yields 

AT-l 

X = -- ($_(o) + $+(0)) ^ (z;- + l;>+) + o{\-^h\ 

The last step is to recall that the sum J2j (-^j + ^j ^) negligible, and thus, 
X = 0{X~^h^). We finally conclude that the error in evaluating the Evans 
function with the Gauss-Legendre method is 

ED = 0{X-'/^h^X-'h\X-''/^h^). 
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